R-matrix analysis of reactions in the ^B compound system 
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Recent activity in solving the 'lithium problem' in big bang nucleosynthesis has focused on the role 
that putative resonances may play in resonance-enhanced destruction of ^Li. Particular attention 
has been paid to the reactions involving the ^B compound nuclear system, d+^Be— >-^B. These 
reactions are analyzed via the multichannel, two-body unitary _R-matrix method using code {EDA) 
developed by Hale and collaborators. We employ much of the known elastic and reaction data, in 
a four-channel treatment. The data include elastic ^He-|-®Li differential cross sections from 0.7 to 
2.0 MeV, integrated reaction cross sections for energies from 0.7 to 5.0 MeV for ^Li(''He,p)*Be* and 
from 0.4 to 5.0 MeV for the ''Li(''He,d)^Be reaction. Capture data have been added to an earlier 
analysis with integrated cross section measurements from 0.7 to 0.825 MeV for ®Li("^He,7)^B. The 
resulting resonance parameters are compared with tabulated values, and previously unidentified 
resonances are noted. Our results show that there are no near d-|-''Be threshold resonances with 
widths that are lO's of keV and reduce the likelihood that a resonance-enhanced mass-7 destruction 
mechanism, as suggested in recently published work, can explain the ^Li problem. 
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I. INTRODUCTION 

Calculations of the abundance of ^Li[l] overestimate 
the value extracted from observations of low-metallicity 
halo dwarf stars 2 , where the stellar dynamics are sup- 
posed to be sufficiently understood to isolate the pri- 
mordial ^Li component. The discrepancy with this (and 
another|31) observation by a factor of 2.2 — >■ 4.2 corre- 
sponds to a deviation of 4.5cr — >■ 5.5(t, a result that has 
only become more severe with time. It is essential to 
determine the nature of this discrepancy as big-bang nu- 
cleosynthesis (BBN) probes conditions of the very early 
universe and our understanding of physical laws relevant 
in an extreme environment. 

Recent attention has focused on the role of reactions 
that destroy A = 7 nuclei at early times < 1 s in the 
big-bang environment [4l |5j . The authors of Ref . 4 , cit- 
ing the TUNL-Nuclear Data Group (NDG) evaluation 
tables[B], (See Table I.) conjecture that the putative | 
resonance near 16.7 MeV may enhance the destruction of 
^Be through reactions like ^Be(d,p)aa and ^Be(d,7)^B if 
the resonance parameters are within given ranges. These 
studies employ the Wigner limit'lO] to determine an up- 
per bound on the contribution of resonances, particularly 
^B, to a resonant enhancement in reactions that destroy 
mass-7 nuclides, ^Li in particular. Because there is a 
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TABLE I. The TUNL-NDG/ENSDF resonances in the ^B 
compound nuclear system [H] for resonances that are low- lying 
with respect to the d-^Be threshold, which occurs at 16.4901 

MeV. 



paucity of data in the region near the d±'^Be threshold 
where the assumed | ^B resonance inhabits, we won- 
dered if the existing data may indicate the presence of 
such a resonance if a multichannel, unitary i?-matrix eval- 
uation is pursued. 

Our motivation for the present study of the ^B com- 
pound system is two-fold. The continuing light nuclear 
reaction program at Los Alamos National Laboratory, T- 
2 Theoretical Division provides light nuclear data for an 
array of end users, including the ENDF and ENSDF com- 
munities. Moreover, we are interested in updating the 
evaluation of the ^B compound system to address the key 
question outlined above for BBN: does a resonance near 
the d+^Be threshold cause an increase in the destruc- 
tion of mass-7 nuclides in the early universe and possibly 
explain the ^Li overprediciton problem? 
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II. THE i?-MATRIX FORMALISM AND EDA 
CODE 

The i?-matrix approach [7H5] is a unitary, multichan- 
nel parametrization that has proven useful for an array 
of nuclear reaction phenomenology, particularly for light 
nucleiTTj. We give only a brief description here and refer 
to the literature for a more complete description |T^ [T5] . 

We consider only 2 — ?> 2 body scattering and reaction 
processes for light nuclear systems. Configuration space 
is partitioned into an interior, strongly interacting region 
and an exterior. Coulomb or non-polarizing interaction 
region by giving a channel radius, ac for each two-body 
channel. The boundary of separation of these regions is 
the channel surface, S = J2c^c- 

The R matrix is computed as the projection on channel 
surface functions, |c) of the Green's function, Gb = {H + 
Cb-E)-^ 



Rc'c = {c!\H + Cb-E)-^\c) 
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(c1A)(A|c) 
E\ — E 



(1) 



where Cb is the Bloch operator, which accounts for the 
presence of a boundary condition, B on the channel 
surface. The Bloch operator ensures that the operator 
H + Cb is a compact, Hermitian operator having a real, 
discrete spectrum. The i?-matrix parameters. Ex and 
l\c — (c|A) describe the spectrum and residues of the 
resolvent operator; they are treated as parameters ad- 
justed to fit the observed data. Both hadronic and elec- 
tromagnetic (ie. 7-(-^B) channels can be handled in this 
approach. The transition matrix, T the square of which 
gives the observables (cross section, etc.) of the theory, 
is given as 



T = pi/^O-iRiO'V'/' - FO-i, 



(2) 



where Rl == (R^^ - L + B)"\ L = pO'O^, and 
F — Im O, where O is the diagonal matrix of outgoing 
(Coulomb) wave functions in the exterior region. 

The i?-matrix approach is implemented in the EDA 
(Energy Dependent Analysis) code developed by Hale 
and collaborators JJJ. The available two-body scatter- 
ing and reaction data is described by minimization of the 
X^ function with respect to variation of the i?-matrix pa- 
rameters Ex and ^xc- 



III. ANALYSIS AND RESULTS 

The i?-matrix configuration, constructed for input into 
the EDA code, is given in terms of the included channel 
partitions (pairs), the LS terms for each partition, and 
the channel radii and boundary conditions. Be for each 
channel. 

We have included in the analysis the hadronic chan- 
nels: d-|-^Be partition with threshold of 16.5 MeV with 
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FIG. 1. Elastic scattering data from [T^ plotted against the 
7?-matrix fit (solid curve) for center-of-niass differential cross 
section vs. '^He lab energy. 



up to Z?-waves, '^He-f ^Li at 16.6 MeV up to P-waves, and 
p+^Be* at 16.7 MeV up to P-waves. The channel radii 
were constrained to lie in the range between 5.5 fm and 
7.5 fm for these. The electromagnetic 7-l-^B channels in- 
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eluded were E^ , j"i , j"i 
with a channel radius of 50.0 fm. 

The ^B analysis is based upon data gathered from 
the literature and stored in the EXFOR/CSISRS 
database[T3]. We include elastic differential cross section 
data for the ^Li-l-'^He channel given in the range of '^He 
lab energy, 1.30 MeV<S(3He)<1.97 MeV[l5]; integrated 
cross section data for ^Li(^He,p)*Be*[T2] where the final 
state channel is an average of the excited-states of the 
quasi-two-body final state of p-|-*Be* given, in the range 
0.66 MeV < E{^YLe)< 5.00 MeV; integrated cross sec- 
tion for the ^Li(3He,d)'^Be[n] in the range 0.42 MeV < 
£;(^He)< 4.94 MeV; and capture data from the "^Li-l-^He 
initial state in the energy range 0.7 MeV<i?(^He)< 0.825 
MeV[IHj. 

Using about 40 parameters, the results of the x^ mini- 
mization result in a r matrix which gives the solid curves 
appearing in Figs. 1-4, plotted along with the data ob- 
tained from references cited in the paragraphs above. The 
fit quality is fair with x^/datum of 1.91, 0.55, 2.38, and 
0.37 for Figs, 1-4, respectively. The fit to the capture 
data, Fig, 4 has been folded with a Gaussian acceptance 
function whose width is 5 keV to match the quoted energy 
resolution in Ref.|18). 

The present i?-matrix parametrization gives a reso- 
nance structure as presented in Table II, The resonance 
poles of the T matrix are determined by diagonalization 
of the complex "energy-level" matrix 
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E^{MeY) r r(keV) Re£o(MeV) Im£o(keV) £(='He)(MeV) Strength 



16.4754 


1/2- 


768.46 


-.1369 


-384.2 


-0.2054 


0.06 


weak 


17.1132 


1/2- 


0.14 


0.5109 


-0.07 


0.7664 


1.00 


strong 


17.2012 


5/2- 


871.63 


0.5989 


-435.8 


0.8984 


0.40 


weak 


17.2809 


3/2- 


147.78 


0.6785 


-73.9 


1.0178 


0.77 


strong 


17.6754 


5/2+ 


33.33 


1.0631 


-16.7 


1.5947 


0.98 


strong 


17.8462 


7/2+ 


2036.21 


1.2339 


-1018.1 


1.8509 


0.15 


weak 


17.8577 


3/2- 


42.52 


1.2454 


-21.3 


1.8681 


0.97 


strong 


18.0582 


3/2+ 


767.11 


1.4459 


-383.6 


2.1689 


0.54 


weak 


18.4229 


1/2+ 


5446.32 


1.8206 


-2723.2 


2.7309 


0.03 


weak 


18.6872 


1/2- 


10278.41 


2.0749 


-5139.2 


3.1124 


0.15 


weak 


19.6192 


3/2- 


1478.22 


3.0069 


-739.1 


4.5104 


0.52 


weak 



TABLE 11. The resonance structure determined in the present 4-channel fit to data as described in the text. The table displays 
the pole location along with J'^ and pole-strength information, as described in the text. 
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Integrated cross section 
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FIG. 2. Reaction data from il6| . integrated cross section 
^He lab energy. 



_ _ , O is the outgoing- 
is the boundary condi- 
Or.. Details are given in 



where Lc{E) = r^idO /drc)0-^[ 
Coulomb wave function, and Be 
tion given at the channel radius, 
Ref.[19j. 

The first column of Table II gives the real part of the 
pole position, Eq = E^ — iT/2, where Eq is one of the 
eigenvalues of the energy-level matrix, Eq. ( 3 ) relative the 
ground state of ^B. The spin-parity is given in the second 
column. The width F is the center-of-mass width in keV. 
The following column restates the real part of the reso- 
nance pole position relative the '^He-l-^Li threshold in the 
center-of-mass. The column labeled i?(^He) is the corre- 
sponding lab energy. The 'Strength' function is the ratio 
of the sum of the channel widths (defined in Ref. [19] ) di- 
vided by the total width, F-^ J2c ^c- Resonances labeled 
'strong' are clearly seen in at least one of the figures. 

The resonance structure shown in Table II differs sig- 
nificantly from that in Table I Possible reasons for the 
discrepancy include the fact that the current analysis is 
the first, to our knowledge, that includes much of the 



available data in the region below £^('^He)< 3.0 MeV in a 
two-body unitary analysis. Several deductions about the 
resonance structure in the TUNL/ENSDF tables rely on 
associated production of ^B experiments and single-level 
i?-matrix parametrizations^. While more data, partic- 
ularly polarization observables, would constrain the cur- 
rent fit with greater confidence, the present analysis ap- 
pears to be the most comprehensive available that ac- 
counts for the available data in a two-body unitary way. 
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FIG. 3. Reaction data from [T7], integrated cross section 
''He lab energy. 



Returning the problem of the overprediction of ^Li in 
current treatments of BBN, we see that the requirements 
for a near-threshold resonance of Refs.[4! and I5i] are dif- 
ficult to arrange given the resonance structure of Table 
IF Both of these works require a narrow resonance, a few 
lO's of kcV in width with a 100 kcV of the ^He-h^Li (that 
is, 200 keV within the d-|-'^Be) threshold in order to ex- 
plain the overproduction of ''Li in BBN reaction network 
codes [1]. 

The current study does not conclusively eliminate the 
possibility of the mechanism of resonant enhancement 
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''Li(%e,Y)''B 
Integrated cross section 
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FIG. 4. Capture data from [TS], integrated cross section vs. 
^He lab energy. 



of mass-7 destruction. The ^B compound system was 
identified originally by Cyburt and Pospelov[l] as play- 
ing a potential role in the destruction of ^Be precisely 
because there isn't much data in the region near the 
d+^Bethreshold. Our analysis is performed on essentially 
the same data that the existent TUNL-NDG analyses[S] 
a were performed, with the smallest energy probed about 
400-500 keV above the ^He+^Li threshold. It might, 
therefore, be suspected that the present data set would 
give no indication of such a low-lying resonance. Our ex- 
perience with R matrix analysis indicates, however, that 
a resonance of lO's of keV in width would likely - but not 
certainly - have contributions 'in the tail' to the observ- 
ables considered in the present study, particularly in the 



^Li('^He,d)^Be integrated cross section of Fig. 3 



IV. SUMMARY, FINDINGS AND FUTURE 
WORK 



We've studied a possible resonant enhancement of the 
destruction of mass-7 (^Be, in particular) in BBN sce- 
narios. The near threshold, narrow state anticipated 
in Refs.fT and j5] appear not to be supported by our 
multichannel, two-body unitary i?-matrix analysis. We 
have reviewed the i?-matrix method implemented in 
the Los Alamos reaction code for light nuclei, EDA 
and have discussed the included data from four chan- 
nels: elastic ^He-t-^Li, 6Li(3He,p)8Be*, 6Li(3He,d)^Be 
and ^U{^Re,-/fB. 

Our analysis determines a resonance structure sig- 
nificantly different from that published in the TUNL- 
NDG/ENSDF compilation 6J , as can be seen by compar- 
ing the results from the present analysis in Table II with 
the table. Table I for the TUNL-NDG/ENSDF analysis. 
Our immediate objective is to incorporate the *Be* fi- 
nal states for each excited state (rather than averaging 
their contribution as we have done in the present analy- 
sis) . This will allow the extension of the present analysis 
to higher energies and the incorporation of polarization 
dataj20j 21, that we have neglected. 

Our findings for the role of a putative resonance in ^B 
near the d-|-^Be threshold as envisioned in Refs.[4l and |5] 
is that their particular mechanism of resonant enhance- 
ment of mass-7 destruction is an unlikely explanation to 
the ^Li problem in BBN, though low-energy data would 
allow a more conclusive statement of this finding or its 
converse. 

This work was carried out under the auspices of the 
National Nuclear Security Administration. 
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